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Abstract:  Tungsten-based  materials  have  been  proposed  as  replacements  for  depleted 
uranium  in  armor-penetrating  munitions  and  for  lead  in  small-arms  ammunition.  A  recent 
report  demonstrated  that  a  military-grade  composition  of  tungsten,  nickel,  and  cobalt 
induced  a  highly-aggressive,  metastatic  rhabdomyosarcoma  when  implanted  into  the  leg 
muscle  of  laboratory  rats  to  simulate  a  shrapnel  wound.  The  early  genetic  changes  occurring 
in  response  to  embedded  metal  fragments  are  not  known.  In  this  study,  we  utilized  two 
cultured  rodent  myoblast  cell  lines,  exposed  to  soluble  tungsten  alloys  and  the  individual 
metals  comprising  the  alloys,  to  study  the  genotoxic  effects.  By  profiling  cell  transcriptomes 
using  microarray,  we  found  slight,  yet  distinct  and  unique,  gene  expression  changes  in  rat 
myoblast  cells  after  24  h  metal  exposure,  and  several  genes  were  identified  that  correlate 
with  impending  adverse  consequences  of  ongoing  exposure  to  weapons-grade  tungsten 
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alloy.  These  changes  were  not  as  apparent  in  the  mouse  myoblast  cell  line.  This  indicates  a 
potential  species  difference  in  the  cellular  response  to  tungsten  alloy,  a  hypothesis  supported 
by  current  findings  with  in  vivo  model  systems.  Studies  examining  genotoxic-associated 
gene  expression  changes  in  cells  from  longer  exposure  times  are  warranted. 

Keywords:  genotoxic;  tungsten;  military;  fragments;  shrapnel;  cancer 


1.  Introduction 

Retained  or  embedded  fragment  wounds  (shrapnel)  have  long  been  a  hazard  of  war.  Modem 
warfare  and  its  associated  technological  advances  have  produced  munitions  that  have  increased  in  both 
range  and  lethality.  So  too  have  the  devices  designed  to  protect  the  warfighter  from  the  effects  of  these 
same  munitions.  Where  previously  an  injury  from  exploded  ordnances  would  likely  have  proved  fatal, 
advances  in  both  protective  equipment  and  medical  treatment  have  rendered  even  the  most  severe 
injuries  survivable. 

One  practice,  that  of  removing  retained  or  embedded  fragments,  has  changed  little  over  time. 
Military  medical  practice  has  long  held  that  fragments  posing  no  immediate  threat  to  the  individual  can 
remain  in  place.  However,  the  composition  of  munitions  has  changed,  and  many  modem  munitions 
may  pose  a  health  hazard  via  long-term  exposure  to  the  retained  embedded  fragments. 

What  has  also  changed  is  the  science  and  study  of  genomics  and  molecular  biology,  including 
research  on  cellular  changes  down  to  the  level  of  genetic  composition  and  alteration.  The  body  of 
knowledge  and  evidence  related  to  genetic  changes  and  disease  processes  is  expanding  at  an 
exponential  rate.  Where  previous  science  could  only  postulate  on  correlative  associations,  emerging 
research  in  genomics  is  now  providing  causative  evidence  of  oncogenic  changes  initiated  by  exposure 
to  various  substances.  Additionally,  as  military  and  civilian  personnel  increasingly  face  survivable 
injuries  and  embedded  fragments  composed  of  potentially  harmful  materials,  it  is  critical  that 
healthcare  providers  be  aware  of  current  policy  and  practice. 

Many  metals  occur  naturally  in  the  environment  and  are  found  in  the  human  body  at  sub-toxic 
levels.  However,  an  increase  in  metal  exposure  may  cause  cellular  changes,  activate  oncogenes  or 
inactivate  tumor  suppressor  genes.  Common  industrial  exposure  to  heavy  metals  can  occur  by 
ingestion,  inhalation  and  dermal  absorption.  However,  in  a  combat  setting,  exposure  can  also  occur  via 
tissue  penetration  by  metal  fragments.  Although  many  metals  are  found  in  combat  situations  that  could 
potentially  result  in  embedded  fragment  wounds,  two  recently  used  materials  are  of  major  concern  to 
the  U.S.  Department  of  Defense. 

Depleted  uranium  (DU)  was  once  the  metal  of  choice  for  U.S.  military  kinetic-energy  munitions. 
DU  saw  its  first  widespread  combat  use  in  the  First  Persian  Gulf  War.  A  number  of  “friendly-fire” 
incidents  left  several  U.S.  personnel  with  wounds  containing  embedded  DU  fragments.  As  a  result  of 
the  unique  chemical  and  radiological  properties  of  DU,  concern  was  raised  over  the  long-term  health 
effects  of  these  embedded  fragments.  A  1993  study  conducted  by  the  Armed  Forces  Radiobiology 
Research  Institute  at  the  request  of  the  U.S.  Army  Surgeon  General  reviewed  the  potential  hazards  of 
embedded  DU  and  concluded  there  were  “sufficient  uncertainties  about  the  long-term  chemical  and 
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radiological  effects  of  embedded  DU  to  warrant  the  initiation  of  animal  toxicology  pilot  studies”  [1]. 
Since  that  time,  a  large  body  of  research  has  evaluated  the  potential  long-term  health  effects  of 
embedded  DU  fragments  [2-20], 

Nonetheless,  concern  still  exists  over  the  long-term  health  and  environmental  effects  of  depleted 
uranium.  As  a  result,  a  search  for  a  suitable  replacement  for  DU  in  armor-penetrating  munitions  was 
undertaken  and  several  tungsten-based  materials  were  identified.  However,  this  was  not  the  first 
combat  use  of  tungsten.  As  early  as  the  1950s,  the  U.S.  military  had  used  tungsten  carbide  in 
armor-piercing  munitions.  As  armor  improved,  tungsten  carbide  was  replaced  with  an  alloy  of  97.5% 
tungsten  and  2.5%  binder  material  with  a  density  of  18.5  g/mL  [21].  Today,  the  composition  of 
weapons-grade  tungsten  alloy  varies,  but  some  of  the  most  widely  used  formulations  are  91%— 93% 
tungsten,  3%-5%  nickel,  and  either  2%-4%  cobalt  or  iron  [22], 

For  years,  tungsten  was  considered  “relatively  safe”  from  a  toxicological  perspective  [23].  However, 
using  a  rodent  model  system  designed  to  investigate  the  health  effects  of  embedded  metal  fragments, 
Kalinich  et  al.  demonstrated  that  a  military-grade  tungsten  alloy  comprised  of  tungsten,  nickel,  and 
cobalt  induced  aggressive  metastatic  rhabdomyosarcomas  in  the  leg  muscles  of  F344  laboratory 
rats  [24],  A  study  characterizing  the  aerosols  produced  in  a  vehicle  struck  by  a  tungsten  penetrator  also 
raised  concerns  over  the  health  effects  associated  with  exposure  to  the  aerosolized  metal  particles  [25]. 
The  authors  concluded  that  further  research  was  needed  to  determine  the  extent  of  risk. 

Embedded  fragments  are  generally  not  surgically  removed  unless  they  pose  a  future  health  risk. 
However,  the  findings  of  Kalinich  et  al.  [24]  raises  new  concerns  about  long-term  health  effects  of 
retained  tungsten  alloy  fragments.  As  a  result,  the  U.S.  Department  of  Defense  (DOD)  mandated  that 
all  fragments  removed  from  DOD  personnel  be  analyzed  for  metal  content  [26],  However,  questions 
still  remain  on  which  embedded  fragments  pose  a  health  risk  and  which  can  remain  in  place  without 
concern.  This  study  seeks  to  address  some  of  those  concerns  with  respect  to  embedded  tungsten  alloy 
fragments.  The  research  described  herein  is  designed  to  explore  and  assess  the  genotoxic  effects  of 
weapons-grade  tungsten  alloy  on  living  cells.  The  purpose  is  to  identify  early  genetic  changes  induced 
by  exposure  to  weapons-grade  tungsten  alloy  that  may  be  used  as  an  indicator  of  future  carcinogenicity. 
Such  information  could  provide  healthcare  providers  with  the  knowledge  required  to  make  informed 
treatment  decisions  when  dealing  with  embedded  fragment  injuries.  In  addition,  this  information  will 
guide  any  changes  in  policy  or  clinical  practice  on  the  removal  of  embedded  fragments  in  general  and 
weapons-grade  tungsten  alloy  fragments  in  particular.  In  this  in  vitro  study,  we  utilized  cultured  rodent 
muscle  cell  lines  to  assess  the  potential  genotoxic-associated  gene  expression  changes  of  the  metals 
comprising  two  military-relevant  tungsten  alloys. 

2.  Experimental  Section 

2.1.  Cell  Culture 

The  rat  myoblast  cell  line  L6  (CRL  #1458)  and  mouse  myoblast  cell  line  C2C12  (CRL  #1772)  were 
both  purchased  from  the  American  Type  Culture  Collection  (Manassas,  VA,  USA).  Cells  were  grown 
as  a  monolayer  in  Dulbecco’s  Modified  Eagle’s  Medium  (DMEM,  Invitrogen,  Grand  Island,  NY,  USA) 
supplemented  with  10%  fetal  bovine  serum  (Invitrogen)  and  1%  penicillin-streptomycin  (Invitrogen)  at 
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37  °C  in  a  humidified  atmosphere  of  5%  CO2  in  air.  Cultures  were  routinely  fed  every  3-4  days,  and 
passaged  weekly.  Experiments  were  conducted  using  cells  at  80%-90%  confluency  as  determined  by 
microscopic  observation. 

2.2.  Test  Metals  and  Exposure  Groups 

Soluble  test  metals  included  the  following  compounds:  C0CI2  6H2O,  NiCf'bFfO,  FeCl3-6H20, 
Na2W04-2H20,  and  TaCls.  All  compounds  were  purchased  from  Sigma-Aldrich  (St.  Louis,  MO,  USA). 
Because  the  WNiCo  and  WNiFe  compounds  are  not  commercially  available,  they  were  reconstituted 
in  the  laboratory  at  91%  W,  6%  Ni,  and  3%  Co  or  Fe  using  the  method  of  Miller  et  al.  [22],  All  metal 
solutions  were  prepared  in  sterile  DMEM  immediately  prior  to  use.  After  thorough  mixing,  the  metal 
solutions  were  added  to  the  cells  to  a  final  metal  concentration  of  10  pg/mL.  There  were  eight 
treatment  groups:  control  (no  metal),  Ta,  W,  Ni,  Co,  Fe,  WNiCo,  and  WNiFe.  Cells  were  treated  for  24  h 
under  the  conditions  described  above.  Metal  treatments,  under  these  experimental  conditions,  did  not 
adversely  affect  the  cells  as  determined  by  microscopic  observation  and  previously  published  studies  [27]. 

2.3.  RNA  Isolation 

After  metal  exposure  for  the  indicated  duration,  the  medium  was  removed  and  the  cells  were 
washed  twice  with  Dulbecco’s  Phosphate  Buffered  Saline  (Invitrogen).  Total  RNA  was  isolated  from 
the  cells  using  commercially  available  kits  (RNeasy  Mini  Kit;  Qiagen,  Valencia,  CA,  USA)  and 
following  the  manufacturer’s  protocol.  Total  RNA  concentration  and  purity  were  determined  by 
measuring  the  absorbance  at  260  and  280  nm  using  a  NanoDrop  ND-1000  spectrophotometer 
(NanoDrop  Technologies,  Inc.,  Wilmington,  DE,  USA).  RNA  samples  with  a  260/280  and  260/230 
ratio  of  greater  than  1.9  were  used  for  microarray  analysis.  Samples  were  stored  at  -80  °C  until  needed. 

2.4.  Microarray  Analysis 

RNA  samples  were  sent  to  the  Genomics  Core  of  the  Lemer  Research  Institute  at  the  Cleveland 
Clinic  (Cleveland,  OH,  USA)  for  whole  genome  RNA  expression  analysis.  The  integrity  of  the  RNA 
was  assessed  using  the  Agilent  2100  Bioanalyzer  (Santa  Clara,  CA,  USA),  and  samples  with  RNA 
Integrity  Number  (RIN)  greater  than  9.0  were  then  processed  and  labeled  before  hybridizing  to  the 
RatRef-12  vl.O  Expression  BeadChip  (Illumina,  San  Diego,  CA,  USA)  for  the  rat  RNA  and  the 
MouseRef-8  v2.0  Expression  BeadChip  for  the  mouse  RNA  (Illumina).  BeadChips  were  imaged  using 
a  BeadArray  Reader  (Illumina),  and  raw  data  output  was  generated  using  the  GenomeStudio  software 
package  (Illumina).  Relevant  output  files  were  accessed  via  FTP  from  the  Genomics  Core  file  servers. 
The  dataset  for  this  study  was  generated  using  raw  data  from  the  probe-level,  quantile-normalized 
BeadArray  data  file.  The  dataset  was  analyzed  at  the  level  of  the  individual  probes,  since  multiple 
probes  interrogate  transcripts  from  individual  genes.  This  data  file  was  used,  as  opposed  to  gene  level 
or  non-quantile  normalized  data  files,  because  normalized  data  accounts  for  possible  variations  of 
sample  input  on  the  BeadArray  and  also  allows  for  the  detection  of  various  isoforms  of  genes.  The 
complete  probe  datasets  were  then  evaluated  using  hierarchical  clustering  and  comparative  marker 
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selection  (CMS)  analyses  with  GenePattern  software  [28]  (Broad  Institute,  Cambridge,  MA,  USA)  to 
determine  if  there  were  any  significant  differences  in  expression  between  the  experimental  groups. 

2.5.  Quantitative  Real-Time  PCR 

To  confirm  the  expression  validity  of  the  selected  genes  found  in  the  microarray  data,  RT-PCR  and 
quantitative  PCR  were  run  on  the  selected  genes  and  a  housekeeping  gene  (HPRT).  Primer  design  was 
based  on  the  transcript  RefSeq  number,  and  performed  using  Primer-BLAST  (NCBI,  Bethesda,  MD, 
USA).  Primers  were  designed  to  span  an  exon-exon  boundary  to  prevent  amplication  of  possible 
genomic  DNA  contamination.  The  primers  were  obtained  from  Integrated  DNA  Technologies 
(Coralville,  IA,  USA),  and  maintained  at  4  °C.  RT-PCR  reactions  were  performed  to  verify  the 
production  of  a  single  amplicon  of  specific  size.  The  primers  used  are  listed  in  Table  1. 

Table  1.  Primer  sequences  used  in  qPCR. 


Gene  Symbol  Accession  # 


AP2A1 

NM  001 10751 1.1 

CHRND 

NM_0 19298.1 

ERG1 

NM_0 1255 1.2 

FBXW5 

NM  00 1025730.1 

FN3K 

NM  00 1109051.1 

FOS 

NM  022 197.2 

GALNT2 

NM  00 1106196.1 

HPRT1 

NM  012583.2 

NMT2 

NM  207590.1 

POLR3E 

NM  001108503.1 

Forward  Primer  Sequence 

GCTGTGTCTCGCCTAAGCCGG 
AACCGCAGTTACCCCATTGAGTGG 
AGCGAACAACCCTACGAGCACCT 
GCGAGCACACCGTGCCTACA 
CCGAGAGAC  AC  A  AG  AACT  GT  GGT  C  A 
CCCACGGTGACAGCCATCTCC 
ATCAGGTTGCAGGGCTGCCG 
CGA  AGT  GTT  GGAT  AC  AGGCC  AGACT 
TGGTTCCTCCCCAGGGAGCAC 
GAGAGAGGGAGGCGGCCAATG 


Reverse  Primer  Sequence 

AGAGCCAGGGTGCAGGAACGA 
C  ATT  GGG  AC  ACTGGGGTCC  ACG 
TCTCCACCAGCGCCTTCTCGTTA 
CCGATCTTCAGCCCCACTGGC 
ATGGAGCAGGGCAGGGACAATC 
GC  ACC  AGC  C  ACT  GC  AGGT  CTG 
GTTGCTCCCCACGTGTCGCA 
TTGGCTTTTCCACTTTCGCTGATG 
CGGGGTGGTGCATAACCGTGG 
GAGCCTGCTCTGACTCTGGCCT 


cDNA  was  generated  via  reverse  transcription  (RT)  with  1  pg  of  total  RNA  and  an  iScript 
Reverse  Transcription  kit  (BioRad,  Hercules,  CA,  USA)  using  the  Gene  Amp  PCR  System  9700 
(Applied  Biosystems,  Carlsbad,  CA,  USA).  The  protocol  for  running  the  RT  is  as  follows: 

Step  1:  5:00  at  25  °C; 

Step  2:  30:00  at  42  °C; 

Step  3:  5:00  at  85  °C; 

Step  4:  Hold  at  4  °C. 

The  selected  gene  expression  values  from  the  microarray  analysis  were  validated  by  real-time  qPCR 
using  SsoFast  EvaGreen  Supermix  (Bio-Rad)  and  the  Cl 000  Thermal  Cycler  (BioRad)  with  the 
following  protocol: 

Step  1:  Run  at  98  °C  for  2:00; 

Step  2:  Run  at  98  °C  for  0:05; 

Step  3:  Run  at  55  °C  for  0:02; 

Step  4:  Repeat  (x40)  steps  2  and  3 

Step  5:  Determine  melt  curve  from  65-95  °C  in  0.5  °C  increments  for  0:05  each. 
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The  amount  of  transcript  measured  in  each  sample  was  normalized  to  the  housekeeping  gene 
HPRT1.  Comparative  values  of  the  control  group  were  used  to  calculate  fold  change,  and  results  are 
expressed  as  the  fold  change  of  normalized  level  of  mRNA  in  the  treatment  over  the  value  obtained  in 
normal  controls. 

3.  Results 

Previous  research  showed  that  embedded  pellets  of  a  military-grade  alloy  of  tungsten/nickel/cobalt 
induced  highly-aggressive  metastatic  rhabdomyosarcomas  in  laboratory  rats  [25],  while  one  composed 
of  tungsten/nickel/iron  did  not  [29],  To  determine  the  genotoxic-associated  transcriptome  alteration 
potential  of  these  materials,  we  tested  soluble  forms  of  the  alloys,  as  well  as  the  individual  metals, 
using  cultured  rat  and  mouse  myoblast  cell  lines  and  commercially  available  gene  chips.  The  results 
obtained  were  compared  to  untreated  controls  and  tantalum-treated  cells.  Tantalum  has  been  shown  to 
be  an  inert  non-reactive  metal  that  has  been  used  in  medical  implants  [30,31]. 

3.1.  Metal  Treatment  of  Myoblast  Cell  Lines 

Prior  to  initiation  of  this  study,  the  effect  of  metal  concentration  and  exposure  time  on  the  L6  and 
C2C12  cells  was  determined.  As  previously  shown  by  others  [27],  metal  concentrations  of  less  than 
100  pg  of  metal/mL  of  growth  medium  with  exposure  times  of  72  h  or  less  did  not  affect  the  growth  or 
viability  of  the  L6  and  C2C12  myoblast  lines.  To  avoid  confounding  effects  of  cytotoxicity,  a  metal 
concentration  of  10  pg/mL  and  a  24  h  exposure  time  were  selected  for  use  in  this  study. 

3.2.  Probe  Characterizations 

The  Illumina  RatRef-12  array  contains  22,523  probes  on  the  bead  chip,  representing  21,910  unique 
genes.  BeadArray  technology  allows  for  a  statistical  determination  of  accurate  detection  for  each 
transcript  by  calculating  signal  intensity  variance  for  replicate  beads  [32],  We  hypothesized  that  the 
low  concentration  and  short  duration  exposures  would  not  result  in  a  robust  change  in  molecular 
phenotype  and  significant  increase  in  the  number  of  genes  expressed  in  the  differentially  treated  cells. 
Therefore,  we  determined  the  absolute  number  of  expressed  genes  in  each  exposure  condition.  In  this 
study,  8,193  probes  were  determined  to  be  accurately  detected  in  all  eight  exposure  groups  at  p<  0.05 
(36%  of  all  probes).  A  detection  /?-value  of  <0.05  signifies  a  5%  possibility  that  the  resulting  signal 
value  for  a  given  feature  (probe)  is  inaccurate.  The  8,193  probes  corresponded  to  8038  unique  genes 
accurately  detected  at  p<  0.05.  By  comparison,  raising  the  p- value  cut  off  to  <0.01  (a  1%  possibility  of 
inaccurate  detection)  resulted  in  7,204  probes  (representing  6,960  unique  genes)  being  accurately 
detected  in  all  eight  exposure  groups.  In  order  to  determine  whether  or  not  metal  exposure  dramatically 
changed  the  transcriptome  (RNA)  of  the  cells,  the  total  number  of  differentially  expressed  rat  genes 
was  analyzed,  and  no  significant  differences  were  observed  in  a  comparison  of  the  number  of 
differentially  expressed  genes  in  each  of  the  eight  sample  groups  between  the  /7-value  cutoff  of  <0.05 
or  <0.01  (Figure  1  and  Table  2).  This  suggests  only  a  small  effect  in  all  treatment  groups  from  metal 
exposure  after  24  h  and  early  gene  expression  changes  are  occurring  while  molecular  and  cellular 
phenotypes  are  not  robustly  altered. 
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Figure  1.  Rat  data  comparison  between  actual  numbers  of  features  accurately  expressed  in 
each  exposure  group  at  the  detection  p-value  <  0.05  and  <0.01. 
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Table  2.  Rat  data  comparison  between  actual  numbers  of  features  accurately  expressed  in 
each  exposure  group  at  the  detection  p-value  <  0.05  and  <0.01. 


Cx 

Ta 

Fe 

Co 

W 

Ni 

WNiCo 

WNiFe 

p  <  0.05 

9482 

9652 

9507 

9647 

9436 

9640 

9526 

9507 

p<0.01 

7701 

8098 

7629 

7750 

7799 

7806 

8020 

7829 

The  Illumina  MouseRef-8  (v2.0)  array  contains  25,697  probes  representing  over  19,100  unique 
genes.  On  these  10,048  probes  (39%  of  all  probes)  were  accurately  detected  in  all  eight  exposure 
groups  at  a  detection  p  value  of  <0.05,  which  corresponded  to  7578  unique  genes.  At  p  <  0.01,  9048 
probes  were  accurately  detected,  corresponding  to  6902  unique  genes.  No  significant  differences  were 
observed  in  a  comparison  of  the  number  of  differentially  expressed  genes  between  the  eight  sample 
groups  (Figure  2  and  Table  3). 


Figure  2.  Mouse  data  comparison  between  actual  numbers  of  features  accurately  expressed 
in  each  exposure  group  at  the  detection  p-value  <  0.05  and  <0.01. 
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Table  3.  Mouse  data  comparison  between  actual  numbers  of  features  accurately  expressed 
in  each  exposure  group  at  the  detection  yO-value  <  0.05  and  <0.01. 


Cx 

Ta 

Fe 

Co 

W 

Ni 

WNiCo 

WNiFe 

p  <  0.05  11,643 

p<0.01  10,103 

11,483 

9695 

11,469 

9953 

11,608 

10,157 

11,772 

10,142 

11,496 

9990 

11,555 

9843 

11,533 

10,008 

3.3.  Hierarchical  Clustering  Analysis  of  Exposure  Groups 

Complete  transcriptome  signatures  can  be  clustered  to  determine  similarity  in  cellular  response  to 
external  stimuli,  such  as  small  molecular  exposure.  We  hypothesized  that  metal  exposure  with  similar 
tumorigenic  outcomes  would  be  associated  with  similar  transcriptome  signatures  as  profiled  by 
BeadArray.  Thus,  we  used  unbiased  hierarchical  clustering  of  gene  expression  levels  from  accurately 
detected  features  of  the  BeadChip  to  classify  exposure  samples  by  groups.  Hierarchical  clustering 
analysis  of  the  mouse  data  produced  a  dendrogram  showing  groupings  of  like  metals,  with  Ta  and  the 
control  group  on  the  same  branch,  W  and  Ni  one  branch,  Co  and  WNiCo  on  one  branch,  and  Fe  and 
WNiFe  on  a  separate  branch  (Figure  3).  While  the  variable  groupings  appear  consistent  with  known 
effects,  the  short  length  of  the  branches  suggests  the  differences  between  the  groups  are  small  and  that 
there  was  not  a  large  effect  of  metal  exposure  after  24  h.  Ratio  comparisons  done  between  Ta  and 
WNiCo  groups  from  the  mouse  microarray  data  did  not  show  fold  changes  greater  than  1.35  among 
features  accurately  detected  in  all  eight  groups,  further  supporting  the  suggestion  that  there  was  not  a 
significant  difference  between  the  exposure  groups  after  24  h. 


Figure  3.  Dendrogram  from  hierarchal  clustering  analysis  of  mouse  microarray  data. 
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The  results  of  hierarchal  clustering  on  the  rat  dataset  showed  that  after  24  h  of  exposure,  metal 
groups  with  similar  tumor-producing  properties,  such  as  WNiCo  and  Ni,  and  Ta  and  the  untreated 
control  group  (Cx)  aligned  on  near-branches  of  the  dendrogram.  Interestingly,  exposure  groups 
believed  to  have  dissimilar  genotoxic  properties,  such  as  W  and  Fe  aligned  on  near  branches  (Figure  4). 
According  to  the  dendrogram,  these  groupings  suggest  that  W  alone  and  Fe  alone  are  more  different 
than  WNiCo  and  Ni  alone.  It  also  suggests  that  effects  observed  in  the  WNiCo  group  may  be  attributed 
more  to  Ni  than  W  after  24  h  exposure. 

The  rat  dendrogram  (Figure  4)  revealed  two  distinct  group  parings  on  different  arms,  Ta  and  Cx  on 
one  arm,  and  Ni  and  WNiCo  on  another  arm.  Additionally,  the  longer  branch  arms  indicate  a  greater 
effect  of  the  metals  on  the  rat  cells  than  was  observed  in  the  mouse  samples.  The  WNiFe  exposure 
group  was  a  distinct  outlier  from  all  other  groups  in  that  the  gene  expression  of  this  group  was 
distinctively  dissimilar  from  all  other  groups. 
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Figure  4.  Dendrogram  from  hierarchal  clustering  analysis  of  rat  microarray  data. 
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The  arrangement  of  the  dendrogram  clustering  was  consistent  for  both  the  rat  and  mouse  datasets, 
regardless  of  the  cutoff  parameters  used  to  run  the  hierarchical  clustering  analysis:  the  complete 
genome,  a  data  set  of  genes  accurately  expressed  in  at  least  one  group,  accurately  expressed  in  6,  7  or 
8  groups,  or  the  dataset  consisting  of  genes  accurately  expressed  in  all  eight  groups.  This  indicates  that 
the  changes  in  gene  expression  observed  between  the  different  exposure  groups  is  consistent  within  the 
dataset,  regardless  of  the  number  of  accurately  detected  probes  used  in  the  analysis. 

3.4.  Comparative  Marker  Selection  and  Ratio  Analysis 

Unbiased  hierarchical  clustering  led  us  to  hypothesize  that  the  most  significant  differences  would 
observed  when  comparing  non-tumor  (Ta  and  Cx)  to  tumor  (W,  Ni)  groups.  However,  while  the 
dendrogram  of  the  mouse  dataset  showed  promise,  comparative  marker  selection  analysis  failed  to 
show  any  significantly  differentially-expressed  genes  (p  <  0.05;  fold  change  greater  than  1.35) 
between  the  two  classes.  Additionally,  ratios  run  between  the  average  signals  of  Tumor  (W,  Ni)  and 
Non-tumor  (Ta  and  Cx)  groups  failed  to  show  any  fold  change  greater  than  1.35  between  the  two  groups. 
The  lack  of  significant  differences  between  the  average  signals  of  accurately  detected  probes  indicates 
a  low  level  of  altered  gene  expression  between  the  groups  following  24  h  exposure  to  the  metals. 

Led  by  the  same  hypothesis  as  for  the  mouse-specific  analysis  and  based  on  results  from  in  vivo 
studies  that  demonstrated  that  WNiCo  and  Ni  exposure  produced  tumors  while  other  metal  exposures 
did  not  [22,25],  two  classes  of  groups  were  created  within  the  rat  dataset:  Tumor  (W,  Ni,  WNiCo),  and 
Non-tumor  (Ta,  Fe,  Cx).  WNiFe  and  Co,  the  two  sample  groups  least  like  the  others  in  the  hierarchical 
cluster  analysis,  were  excluded.  A  comparative  marker  selection  analysis  was  run  to  determine  the 
number  of  differentially  expressed  genes  between  the  two  classes. 

Comparative  marker  selection  analysis  results  were  filtered  for  fold  changes  of  genes  between 
exposure  groups  of  >1.35  with  a  />value  <  0.05.  Results  showed  upregulation  of  10  genes  in  the 
“Non-tumor”  class,  or  NON  (Cx,  Ta,  Fe),  and  14  genes  in  the  “Tumor”  class,  or  TUM  (W,  Ni, 
WNiCo).  The  fold  change,  probe  ID  number,  and  gene  name  are  shown  in  Table  4. 

Additionally,  the  Ta  and  WNiCo  groups  were  selected  to  represent  each  class  in  pairwise  analysis. 
Genes  that  were  accurately  detected  in  both  Ta  and  WNiCo  that  had  a  ratio  indicating  fold  change 
between  the  two  groups  of  >1.35  or  its  inverse  <0.74  were  selected  to  determine  amount  of  fold  change 
between  the  two  groups  and  plotted  (Figure  5).  Nine  hundred  and  sixteen  (916)  genes  met  these  criteria; 
469  that  showed  a  higher  level  of  gene  expression  in  the  WNiCo  group,  and  447  that  showed  a  higher 
level  of  gene  expression  in  the  Ta  group.  If  the  data  points  make  a  straight  line  going  from  the  origin 
out  to  high  x-  and  y-values,  then  the  variables  are  said  to  have  a  positive  correlation.  In  this  case,  the 
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data  points  that  are  outliers  from  the  trendline  are  of  interest,  because  these  points  indicate  that  the 
gene  expression  values  are  substantially  different,  which  is  what  would  be  expected  if  there  is  a 
significant  difference  in  effect  size  between  WNiCo  and  Ta  exposures. 


Table  4.  Results  of  CMS  analysis  between  Non-tumor  class  (NON)  and  Tumor  class 
(TUM),  showing  those  genes  which  showed  a  >  1.50  fold  change  between  the  two  groups. 


Upregulated  in  Class 

Probe  ID 

Gene  Name 

Fold  Change 

NON 

ILMN  1366318 

RGD 1 306926  PREDICTED 

1.68 

NON 

ILMN_1 376639 

DYX1C1 

1.42 

NON 

ILMN_1 359445 

KCNJ9 

1.86 

NON 

ILMN_1352946 

ZBED3 

1.51 

NON 

ILMN_1 373980 

CDGAP  PREDICTED 

1.36 

NON 

ILMN_1 649958 

LOC294726 

1.36 

NON 

ILMN_1 376505 

TTLL1 

1.64 

NON 

ILMN1361716 

LGALS2 

1.36 

NON 

ILMN_1 375469 

RGD  1 307826  PREDICTED 

1.38 

NON 

ILMN  1354019 

SARDH 

1.44 

TUM 

ILMN  1367175 

RGD  1 562265  PREDICTED 

1.37 

TUM 

ILMN_1363142 

NMT2 

1.49 

TUM 

ILMN_1350240 

RGD131 1358 

1.37 

TUM 

ILMN_1353413 

RGD  1308 1 06  PREDICTED 

1.58 

TUM 

ILMN_1358380 

CNOT3PREDICTED 

1.37 

TUM 

ILMN  1360156 

RGD  156271 6  PREDICTED 

1.47 

TUM 

ILMN_1 364447 

MLANAPREDICTED 

1.51 

TUM 

ILMN_1 368669 

RGD  156421  ^PREDICTED 

1.43 

TUM 

ILMN  1368356 

FOS 

1.70 

TUM 

ILMN  1 369005 

EGR1 

1.47 

TUM 

ILMN_1351892 

RGD1 560234  PREDICTED 

1.49 

TUM 

ILMN_1 362498 

ZFP84PREDICTED 

1.40 

TUM 

ILMN_1362589 

CWF 1 9L 1  PREDICTED 

1.36 

TUM 

ILMN 1376212 

ACE2 

1.52 

Figure  5.  Plot  of  WNiCo  and  Ta  signal  values  for  fold  change  >  1.35  and  <  0.74. 


Int.  J.  Environ.  Res.  Public  Health  2014, 11 


2932 


3.5.  PCR  Analysis 

To  validate  the  microarray  results  of  the  rat  cell  samples,  we  utilized  quantitative  real-time  PCR 
analysis,  a  gold-standard  independent  methodology  with  ability  to  detect  minimal  fold-changes  in 
expression.  Six  genes  were  selected  from  the  microarray  data  that  showed  a  fold  change  greater  than 
1.35  between  Ta  and  WNiCo;  three  with  greater  expression  in  the  Ta  group  (AP2A1,  CHRND, 
FBXW5),  and  three  with  greater  expression  in  the  WNiCo  group  (FN3K,  GALNT2  POLR3E)  (Table  5). 
Additionally,  three  genes  (FOS,  EGR1  and  NMT2)  were  selected  that  were  among  the  CMS  results. 
All  but  two  of  the  genes  (CHRND  and  FN3K)  were  accurately  detected  in  all  eight  exposure  groups, 
with  CHRND  accurately  detected  in  six  groups  (all  except  for  W  and  WNiCo  groups),  and  FN3K, 
which  was  accurately  detected  in  three  groups  (Ta,  Fe  and  WNiCo  groups)  (Table  5).  Paired  /-test 
/^-values  for  the  effect  of  WNiCo  and  Ta  were  found  to  be  consistent  with  the  microarray  data.  Four 
genes  (AP2A1,  CHRND,  FN3K  and  GALNT2)  were  differently  expressed  at  statistically  significant 
levels  between  the  exposure  groups  using  2-tailed  /-tests  (p  =  0.03),  and  directionality  of  expression 
differences  were  validated  in  six  of  the  genes  (AP2A1,  CHRND,  FBXW5,  FN3K,  POLR3E).  However, 
for  two  genes  (GALNT2  and  POLR3E),  PCR  expression  directionality  was  reverse  of  the  microarray 
results  (Figure  6).  A  two-tailed  /-test  between  CHRND  expression  in  WNiCo  and  Ta,  and  between 
WNiCo  and  Cx  showed  significance  {p  <  0.03  and  p  <  0.05,  respectively),  indicating  a  difference  in 
expression  between  the  WNiCo  group  and  both  Ta  and  Cx.  AP2A1,  FN3K  and  GALNT2  showed 
similar  results  for  WNiCo  and  Ta  (p  <  0.03,  p  <  0.04,  and  p  <  0.04  respectively).  Based  on  this 
information,  it  was  concluded  that  the  expression  values  and  fold  changes  observed  in  microarray  data 
correlated  with  PCR  results. 


Table  5.  Genes  selected  for  PCR  validation,  with  rat  microarray  average  signal  values. 


ACCESSION 

SYMBOL 

Cx 

Ta 

Co 

W 

Ni 

WNiCo 

WNiFe 

XM_001078388.1 

AP2A1 

45.70 

57.40 

41.10 

48.00 

43.70 

25.30 

36.40 

NM_0 19298.1 

CHRND 

26.90 

22.80 

19.50 

23.80 

10.50 

7.40 

17.50 

NM_00 1025730.1 

FBXW5 

96.10 

74.30 

60.90 

66.50 

75.40 

43.00 

57.40 

XM_573235.2 

FN3K 

5.80 

15.20 

-2.80 

11.70 

9.30 

26.20 

10.10 

XM_00 105408 1.1 

GALNT2 

38.80 

20.20 

35.10 

31.10 

27.70 

50.10 

31.10 

XM_00 1075970.1 

POLR3E 

35.60 

47.00 

67.10 

39.90 

47.40 

75.50 

63.60 

NM_022 197.1 

FOS 

27.50 

19.10 

18.60 

34.20 

44.90 

47.60 

40.00 

NM_0 1255 1.1 

EGR1 

209.60 

253.10 

296.50 

367.60 

429.80 

415.90 

426.60 

NM_207590.1 

NMT2 

15.30 

20.60 

20.30 

37.20 

22.90 

28.70 

29.80 

3.6.  Genes  of  Interest 

Of  the  genes  that  were  selected  from  the  rat  dataset  for  validation  via  qPCR,  several  are  known  to 
be  associated  with  tumorigenic  pathways,  apoptosis,  or  cellular  repair.  Of  particular  interest  were  the 
genes  CHRND,  EGR1,  FN3K,  FOS  and  NMT2  (Table  6). 
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Figure  6.  PCR  expression  values  for  selected  genes  for  the  Ta  control  and  WNiCo  groups. 
Data  are  the  mean  of  three  measurements.  Error  bars  represent  standard  error. 


Table  6.  Validated  genes  found  to  be  associated  with  genotoxic  pathways. 


Gene 

Name 

Accession  # 

CHRND 

Chrnd  cholinergic  receptor,  nicotinic,  delta 

NM_0 19298.1 

EGR1 

Egrl  early  growth  response  1 

NM_0 1255 1.1 

FN3K 

Fructosamine  3  kinase 

NM_00 1109051.1 

FOS 

Fos  FBJ  murine  osteosarcoma  viral  oncogene  homolog 

NM_022 197.2 

NMT2 

N-myristoyltransferase  2 

NM 207590.1 

Nicotinic  cholinergic  receptor  subunit  delta  (CHRND),  was  accurately  detected  in  six  of  eight 
groups,  and  was  differentially  expressed  to  a  greater  extent  in  the  Ta  group  as  opposed  to  the  WNiCo 
group,  with  a  fold  change  of  3.08.  EGR1  was  accurately  detected  in  all  eight  exposure  groups,  with  an 
upregulated  fold  change  of  1.6  more  in  the  WNiCo  group  than  in  Ta  group.  While  the  PCR  validated 
the  expression,  if  not  the  directionality  in  both  the  Ta  and  WNiCo  groups,  /-test  results  showed  no 
statistically  significant  difference  between  the  expression  values  of  two  groups  in  the  PCR  results. 
Numerically,  the  average  expression  value  of  the  Ta  and  Cx  groups  were  very  close  (0.97  and  0.95, 
respectively);  expression  was  slightly  higher  in  WNiCo  (1.02)  than  both  Ta  and  Cx  values. 

Analysis  of  the  microarray  data  showed  FOS  upregulated  in  the  tumor  class  of  the  CMS  with  a  fold 
change  of  1.7,  and  a  fold  change  2.5  times  higher  in  WNiCo  than  Ta.  There  was  no  statistical 
difference  in  expression  between  WNiCo  and  Ta  in  the  PCR  results,  although  both  were  validated  by 
the  test,  and  the  average  of  the  WNiCo  results  were  higher  than  both  Ta  and  Cx. 

FN3K  was  upregulated  in  the  tumor  class  in  the  CMS  data,  and  in  WNiCo  alone  with  a  fold  change 
of  1.7  times  greater  in  WNiCo  than  in  Ta.  A  /-test  between  the  Ta  and  the  WNiCo  group  PCR  results 
showed  statistically  significant  differences  in  the  expression  of  FN3K  in  the  WNiCo  groups  over  the 
Ta  groups  (p  =  0.03). 
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4.  Discussion 

The  purpose  of  this  study  was  to  look  at  very  early  events  occurring  at  the  genetic  level  in  cultured 
rodent  muscle  cells  after  exposure  to  tungsten  alloys  and  their  individual  component  metals,  with  the 
aim  of  identifying  early  changes  that  might  contribute  to  the  development  of  metal-induced  tumors. 
It  was  found  that  after  24  h  of  exposure  to  the  different  metal  groups,  subtle  changes  in  gene 
expression  between  the  exposure  groups  were  observed  in  the  rat  muscle  cells,  but  were  not  seen  in  the 
mouse  muscle  cells.  Differences  in  the  species’  genetic  response  to  metal  exposure  may  account  for 
this,  and  longer  exposure  intervals  would  be  required  to  see  if  changes  in  expression  become  more 
pronounced  over  time. 

Hierarchical  clustering  analysis  yielded  consistent  results  showing  similarities  between  the  basic 
networks  of  samples.  For  the  mouse  dataset,  while  Ta  and  Cx,  and  W  and  Ni,  paired  together  on 
unique  branches,  WNiCo  and  Co  consistently  paired  on  a  separate  branch,  and  Fe  and  WNiFe  paired 
on  a  branch  independent  from  the  others.  For  the  rat  dataset,  Ta  and  Cx,  Ni  and  WNiCo,  and  W  and  Fe 
consistently  showed  up  paired  together  on  different  arms  of  the  dendrogram  hierarchy,  indicating  Ta 
and  Cx  displayed  similarities,  and  were  sufficiently  dissimilar  to  Ni  and  WNiCo,  which  showed 
similarities  to  each  other.  Co  showed  some  similarity  to  Ni  and  WNiCo,  and  they  were  on  the  same 
branch.  WNiFe  was  an  outlier,  showing  up  on  a  completely  separate  arm  from  all  other  groups. 

As  was  indicated  in  the  results  section,  this  information  initially  looked  promising,  but  further 
analysis  of  the  mouse  data  showed  only  very  small  expression  differences  between  the  groups  as 
evidenced  by  fold  changes  of  less  than  1.35  and  the  small  number  of  accurately  detected  probes, 
indicating  that  after  24  h  exposure,  there  was  not  very  much  activity  stimulated  by  exposure  to  the 
metals.  This  may  be  due  to  a  delayed  response  in  mouse  cells  to  exposure  to  the  metals,  a  higher 
tolerance  than  rat  cells  to  metal  toxicity,  or  a  better  adaption  mechanism  than  rat  cells.  Current 
research  [33]  has  shown  that  mice  implanted  with  WNiCo  pellets  eventually  develop 
rhabdomyosarcomas,  but  the  time  to  tumor  formation  is  far  longer  (12-15  months)  than  that  seen  for 
WNiCo-implanted  rats  (3  months)  [25],  This  suggests  that  mice  may  have  a  protective  mechanism  that 
impedes  the  genetic  damage  leading  to  tumorigenesis  until  later  in  their  lifecycle,  at  which  time  it 
either  shuts  off,  or  is  overwhelmed  by  the  metal  exposure.  Whatever  the  cause,  further  study  utilizing 
both  in  vitro  and  in  vivo  model  systems  is  warranted. 

Since  the  differential  changes  between  control  and  metal-treated  mouse  cells  were  not  statistically 
significant,  genes  for  PCR  validation  were  selected  from  the  rat  microarray  data  set  where  probe  IDs 
were  accurately  expressed  in  all  eight  groups.  The  genes  were  selected  for  signal  strength  rather  than 
gene  function  in  order  to  validate  overall  gene  expression  in  the  microarray  data  using  real  time  qPCR, 
considered  the  gold  standard  for  gene  expression  validation.  This  process  was  part  of  the  discovery 
process,  and  helped  to  validate  the  overall  microarray  data.  None  of  the  genes  showed  a  fold  change 
greater  than  2.50,  but  while  the  differences  in  expression  between  the  different  groups  was  subtle,  it  is 
important  to  note  that  changes  were  observed,  indicating  that  there  are  metal-induced  changes 
occurring  in  the  rat  muscle  cells  as  soon  as  24  h  after  exposure. 

Three  genes  were  validated  in  the  qPCR  showing  both  directionality  and  statistically  significant 
differences — AP2A1,  CHRND  and  FN3K.  AP2A1  (adaptor-related  protein  complex)  is  associated 
with  endocytosis,  protein  transport,  cargo  selection  and  vesicle  formation,  but  does  not  appear  to  be 
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directly  associated  with  tumorigenic  pathways  [34],  FN3K  (fructosamine  3  kinase)  is  involved  in 
cellular  repair,  maintenance  and  defense  against  toxins,  primarily  repairing  protein  damage  caused  by 
glucose,  and  is  responsible  for  protein  deglycation,  and  considered  a  housekeeping  gene  in  cellular 
maintenance  [35], 

Interestingly,  CHRND,  a  member  of  the  nicotinic  acetylcholine  receptor  family,  is  a  recognized 
candidate  for  tumor  associated  rearrangements  and  as  one  of  the  muscle-specific  genes  that  is 
expressed  in  rhabdomyosarcomas  [36,37].  It  is  considered  diagnostic  tumor  marker  [36,38].  However, 
while  the  difference  in  CHRND  expression  between  Ta  and  WNiCo  was  statistically  significant,  the 
expression  value  was  higher  in  the  Ta  group  than  the  WNiCo  group.  One  possible  explanation  for  this 
is  that  the  gene  may  initially  be  suppressed  early  in  the  response  to  WNiCo  exposure.  To  evaluate 
adequately,  it  would  be  necessary  to  follow  CHRND  expression  across  longer  exposure  times  or  to 
utilize  tumor  tissue  from  animal  studies. 

Identification  of  genes  within  known  tumorigenic  or  tumor  suppression  pathways  that  either 
up-regulate  or  down-regulate  in  response  to  metal  exposure  is  crucial  in  understanding  the  impact  of 
long-term  exposure  to  embedded  metals,  including  the  tungsten  alloys.  Correlation  of  early  genetic 
changes  in  an  in  vitro  model  with  known  tumorigenic  outcomes  in  an  in  vivo  model  will  provide  the 
information  needed  to  make  informed  decisions  on  the  proper  treatment  of  embedded  metal  fragment 
wounds.  Of  the  genes  identified  in  the  rat  microarray  data  and  validated  via  qPCR,  several  bear  further 
investigation  as  they  may  indicate  initiation  of  oncogenic  processes  that  correlate  with  metal  exposure. 
They  include  EGR1,  FOS,  FN3K,  and  NMT2. 

EGR1  is  essential  for  mediation  of  the  p5 3 -independent  c-Myc-induced  apoptosis.  It  is  also 
implicated  in  progression  of  certain  types  of  prostate  cancers,  and  possibly  implicated  in  tumor 
suppression  [39].  The  expression  of  EGR1  and  the  slight  difference  between  the  tumor  and  non-tumor 
groups  may  allude  to  the  possibility  of  a  greater  difference  in  expression  between  the  groups.  Even  this 
slight  difference  in  expression  of  the  gene  between  these  groups  after  only  24  h  may  indicate  that 
cellular  damage  and  the  initiation  of  oncogenic  pathways  is  occurring  almost  immediately  after 
exposure. 

FOS  (FBJ  osteosarcoma  oncogene)  is  a  canonical  oncogene,  and  is  part  of  the  AP-1  transcription 
factor,  which  upregulates  transcription  of  genes  involved  in  proliferation,  differentiation  and  defense 
against  invasion  and  cell  damage.  It  is  associated  with  cell  transformation  and  progression  of  cancer. 
Interestingly,  a  recent  study  showed  that  the  disruption  of  certain  oncosupression  transcriptors  have 
been  shown  to  induce  both  FOS  and  EGR1  expression  in  the  early  stages  of  tumor  development  [40], 
FN3K  (FNSK)  is  considered  a  housekeeping  gene  that  is  activated  in  response  to  an  increase  in 
interleukin  ip.  It  involved  in  cellular  repair  in  response  to  glucose-related  protein  damage.  Specifically, 
it  catalyzes  phosphorylation  of  fructosamines  that  are  formed  by  glycation,  the  nonenzymatic  reaction 
of  glucose  [35,41],  Expression  of  this  gene  in  the  WNiCo  exposure  group  may  indicate  early  cellular 
damage  that  initiates  a  cellular  repair  process  soon  after  exposure. 

NMT2  is  also  associated  with  oncogenic  pathways.  It  is  part  of  the  gene  group  associated  with 
myristoylation  (irreversible  protein  modification)  [42],  There  is  some  evidence  that  NMT2  may  play  a 
role  in  the  pathogenesis  of  human  colorectal  cancer.  Interestingly,  the  expression  of  NMT2  is  rapidly 
upregulated  within  minutes  in  response  to  inflammatory  signals,  a  component  of  tumorigenesis  [43]. 
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Furthermore,  NMT2  induction  (activation  by  a  molecule  inactivating  a  repressor  gene)  in  direct 
response  to  carcinogen  exposure  (dioxins)  has  also  been  demonstrated  [44], 

5.  Conclusions 

While  the  recent  interest  over  the  effects  of  tungsten  alloys  has  largely  focused  on  their  effects  on 
muscle  tissue,  it  is  becoming  increasingly  evident  that  these  alloys  have  toxic  effects  on  other 
physiologic  systems,  and  gene  expression  profiles  are  consistent  with  some  of  the  results  observed  in 
this  study.  Recent  research  on  the  effects  of  inhaled  WNiCo  showed  a  dose-dependent  inflammatory 
response  that  resulted  in  the  expression  of  genes  associated  with  apoptosis,  DNA  damage  and  repair, 
and  carcinogenesis  [45],  WNiCo  was  also  recently  shown  to  have  cytotoxic  effects  on  human  brain 
and  kidney  cells  at  a  greater  rate  than  W  alone  [46],  Clearly,  additional  research  is  needed  to  assess  the 
cytotoxic  and  genotoxic  effects  of  tungsten  alloy  exposure  occurring  by  various  routes  of  internalization. 
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